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Abstract. STAR has measured a variety of strange particle species in p + p collisions at y/s = 200 GeV. 
These high statistics data are ideal for comparing to existing leading- and next-to-leading order perturbative 
QCD (pQCD) models. Next-to-leading (NLO) models have been successful in describing inclusive hadron 
production using parameterized fragmentation functions (FF) for quarks and gluons. However, in order 
to describe identified strange particle spectra at NLO, knowledge of flavor separated FF is essential. Such 
FF have recently been parameterized using data by the OPAL experiment and allow for the first time to 
perform NLO calculation for strange baryons. In fact, comparing the STAR A data with these calculations 
allow to put a constraint on the gluon fragmentation function. We show that the Leading-order (LO) event 
generator PYTHIA has to be tuned significantly to reproduce the STAR identified strange particle data. In 
particular, it fails to describe the observed enhancement of baryon-to-meson ratio at intermediate pr (2-6 
GeV/c). In heavy-ion (HI) collisions this observable has been extensively compared with models and shows 
a strong dependency on collision centrality or parton density. In the HI context the observed enhancement 
has been explained by recent approaches in terms of parton coalescense and recombination models. 
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1 Introduction 



Perturbative QCD has proven to be successful in describ- 
ing inclusive hadron production in elementary collisions. 
Within the theory's range of applicability, calculations 
at next-to-leading order (NLO) have produced accurate 
predictions for transverse momentum spectra of inclusive 
hadrons at different energy scales [HE]- With the new 
high statistics proton-proton data at y/s — 200 GeV col- 
lected by STAR, we can now extend the study to identified 
strange hadrons as well as strange resonances. 

The perturbative QCD calculation applies the factor- 
ization ansatz to calculate hadron production and relies on 
three ingredients. The non-perturbative parton distribu- 
tion functions (PDF) are obtained by parameterizations 
of deep inelastic scattering data. They describe quantita- 
tively how the partons share momentum within a nucleus. 
The second part, which is perturbatively calculable, con- 
sists of the parton cross-section amplitude evaluated to LO 
or NLO using Feynman diagrams. The third part consists 
of the non-perturbative Fragmentation functions (FF) ob- 
tained from e + + e~ collider data using quark-tagging al- 2 Data Analysis 



Knichl, Kramer and Potter (KKP) have shown that FF 
are universal between e + + e~ and p + p collisions [3]. 

The theoretical mechanisms of baryon production have 
been difficult to understand and different attempts have 
been made [4]. In the string fragmentation approach the 
production of baryons is intimately related to di-quark 
production from strings. They then combine with a quark 
to produce a baryon. In NLO calculations, baryon produc- 
tion is based on the knowledge of baryon fragmentation 
functions (FF) from e + + e~ collisions. So far the only 
baryon FF which has been accurately measured and pa- 
rameterized is that of the proton [5]. Other groups have 
used a statistical approach to calculate FF [6|. 

In the following section, we compare our p+p data to 
PYTHIA, the most commonly used leading-order Monte 
Carlo event generator for elementary collisions. In par- 
ticular, we study predictions for baryons and the ratios of 
baryons to mesons and see how parameter tunes affect the 
data. We then compare our data with more sophisticated 
NLO calculations. 



gorithms. These parameterized functions are sufficiently 
well known for fragmenting light quarks, but less well 
known for fragmenting gluons and heavy quarks. Recently, 
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The present data were reconstructed using the STAR de- 
tector system which is described in more detail elsewhere 
[7]. The main detector used in this analysis is the Time 
Projection Chamber (TPC) covering the full acceptance 
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in azimuth and a large pseudo-rapidity coverage (| 7/ |< 
1.8). A total of 14 million non-singly diffractive (NSD) 
events were triggered with the STAR beam-beam coun- 
ters (BBC) requiring two coincident charged tracks at 
forward rapidity. Due to the particulary low track mul- 
tiplicity environment in p+p collisions, only 76% of pri- 
mary vertices are found correctly; from the remainder, 
14% are lost and 10% are badly reconstructed as a MC- 
study showed. Of all triggered events, 7 million events 
passed the selection criteria requiring a valid primary ver- 
tex within 50cm along the beam-line from the center of 
the TPC. The strange particles were identified from their 
weak decay to charged daughter particles. The following 
decay channels and the corresponding anti-particles were 
analyzed: Kg -> tt+ + 7r" (b.r. 68.6%), A -> p + 7r~(b.r. 
63.9%) ,S~ -> A + TT-(b.r. 99.9%). Particle identification 
of the daughters was achieved by requiring the dE/dx to 
fall within the 3<7-bands of the theoretical Bethc-Bloch pa- 
rameterizations. Further background in the invariant mass 
was removed by applying topological cuts to the decay 
geometry. Corrections for acceptance and particle recon- 
struction efficiency were obtained, as a function of pr, by 
a Monte-Carlo based method of embedding simulated par- 
ticle decays into real events and comparing the number of 
simulated and reconstructed. 



3 Comparison to PYTHIA 

3.1 Strange particle spectra 

One of the most widely used models for simulating ele- 
mentary collisions is PYTHIA [8|. It is a parton-shower 
based event generator that includes leading order parton 
processes and parton fragmentation based on the Lund 
Model. The parton distributions of the initial state pro- 
tons can be chosen from an array of PDFs (here we use 
CTEQ5M). The model is being actively used and the au- 
thors have recently released a version with completely 
overhauled multiple scattering and shower algorithms (ver- 
sion 6.3) [9]. The PYTHIA version used in this paper is 
6.317. 

The string fragmentation based on the Lund Model 
requires only two parameters to define the shape of the 
fragmentation function and is universal for all light quark 
flavors. Baryons are produced from di-quarks and their 
probability is suppressed with respect to qq pair produc- 
tion. Next-to-leading order processes can be "simulated" 
in PYTHIA by tuning the K-factor (MSTP(33)) or by in- 
creasing the parton shower activity. This will enhance the 
relative probability of hard processes of type quark-gluon 
and thus mock-up the contributions from higher order pro- 
cesses. 

In figure [T] (upper row), we compare PYTHIA calcu- 
lations for strange mesons and baryons to the measured 
STAR data. Whereas the default parameters (blue line) 
agree quite well for the Kg, they clearly underestimate the 
yields at intermediate px for the A and S~ . By increasing 
the K-factor to 3 (red line) wc achieve a reasonable agree- 
ment with the data. In figure [T] (lower row), we compare 



Table 1. STAR dN/dy for various baryons from p+p collisions 
at yfs = 200 GeV (\y\ <0.5) compared to PYTHIA 6.317. 
Pythia Baryon tune is denned as PARJ(1)=0.125 (D=0.1) and 
PARJ(3)=0.5 (D=0.4). 



Particle 


STAR dN/dy 


PYTHIA 


PYTHIA tuned 


proton 


0.11 ± 0.01 


0.096 


0.11 


A (FD) 


0.0385 ± 0.0035 


0.0297 


0.0371 




0.0026 ± 0.0009 


0.0020 


0.0029 



PYTHIA to the strange resonances K* , <p and E* . Again, 
only when applying a higher K-factor does the calculation 
agree with our data. 

In summary, PYTHIA is capable of describing px spec- 
tra for a variety of particles from p + p collisions at RHIC 
energies. However, we have presented evidence that a tune 
of the LO K-factor is necessary in particular for strange 
baryons and resonances. Of course, we have not explored 
all possibilities of parameter "tunes" and there may be 
other, equivalent ways of reproducing the data. 

What are the possible reasons for this discrepancy? 
The "naive" reason, supported by the K-factor tune, is 
that higher order contributions may be significant. How- 
ever it is troubling that the pions do not seem to require 
this tune as shown previously [14] . even though a similar 
study of K-factors for non-identified hadrons found that 
at y/s = 200GeV a value of 3 was needed p~5] . 

Another, perhaps more natural explanation, may be 
related to fragmentation functions for baryons in PYTHIA. 
In the next section we will discuss possible changes to the 
baryon production parameters, ie. the di-quark suppres- 
sion factors, which may help solve this discrepancy. 



3.2 Baryon production 

In string models, baryon production in its simplest form 
is understood via the production of di-quark pairs from 
string-breaking and their recombination with other quarks. 
This process is suppressed with respect to q — q pairs from 
string-breaking resulting in systematically lower baryon 
yields than mesons. The default value for the suppres- 
sion factor is P(qq) = 0.1 x P(q). We have increased this 
value to 0.125 (PYTHIA parameter PARJ(l)). Similarly 
the strange di-quarks are suppressed with a default fac- 
tor P(sq) = 0.4 x P(q), which we have increased to 0.5 
(PYTHIA parameter PARJ(3)). 

In table [1] we show recently measured baryon yields in 
p + p collisions at mid-rapidity. Values for A have been 
corrected for feed-down (FD) from ^"-decays. From the 
values in the table, it is clear that the tuned values for 
PYTHIA are in better agreement with the experimental 
measurements of STAR than the default values. However, 
it must be said that the agreement is confined to low 
Pt and that this tune does not change the shape of the 
PYTHIA spectra to improve the high pt part. 
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Fig. 1. (Top) Minimum-bias pt spectra for Kg, A and S at (| y |< 0.5) from p + p at y/s = 200GeV. (Bottom) K* and 0, and 
X 1 * pt spectra at mid-rapidity. In the left panel, black symbols are K°* and blue symbols are K + * [11U12|[T3] . 



3.3 Baryon to Meson ratios 

Recent heavy-ion data from STAR show a large enhance- 
ment of the baryon to meson ratios at intermediate px, 
which is associated with parton coalescence and recom- 
bination models [21]. A and Kg are ideal candidates for 
comparing baryon to meson production at these momenta 
since they can be cleanly identified via the topological re- 
construction method described at the beginning. 

In Figure [2] we show the measured A/Kg ratio vs px 
measured by STAR, together with 3 different calculations 
by PYTHIA. Open symbols depict S~ feed-down cor- 
rected A yields. Clearly, the default PYTHIA calculation 
lies well below the data. Increasing the LO K-factor does 
not improve the ratio much at low pt, although it does 
describe the ratio at high px- However, using the tuned 
baryon parameters discussed in the previous section im- 
proves the agreement at low px considerably. Thus, we 
need to us a combination of both K- Factor and baryon pa- 
rameter tune to simultaneously describe the spectra and 
the ratios. 

This result triggers the interesting question as to the 
possible energy dependence of this baryon production pa- 
rameter. To investigate this further, we have used data for 
strange mesons and baryons from the UA1 collaboration, 
which measured p + p collisions at s/s — 630 GeV and 
produced the particle ratio presented in figure [3] [5D] . 

The figure clearly shows that the disagreement with 
default PYTHIA for the baryon to meson ratio is not spe- 
cific to our energy scale but also exists at higher energies. 
At y/s = 630 GeV, the difference between PYTHIA and 



• STAR, p+p,^=200 GeV 






~~ O STAR (feeddown corrected) 






PYTHIA6.317 






PYTHIA6.317, K=3 






PYTHIA 6.317, PARJ(1)=0.125, PARJ(3)=0.5 






7 liflMm 






- MV \%K 1. 

o / .-' V J B A 

•0/ j> V /V i\ 

- }K/ ^ 






Preliminary Data 




, , , 1 



P T [GeV/c] 

Fig. 2. Ratio of A+A/2 x Kg vs p T from STAR data compared 
to three different tunes of PYTHIA 6.317. Data are shown with 
and without feed-down correction, whereas PYTHIA calcula- 
tions are corrected for feed-down. 



data is about a factor of 3 and the enhancement of ^/Kg 
is twice as large as in STAR. Even when tuning PYTHlA 
to the same values as for STAR the discrepancy between 
data and model remains large. This may be an indication 
that the effects observed in this ratio in heavy-ion data 
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Fig. 3. Ratio of A/Kg vs pr from UA1 data compared to two 
different tunes of PYTHIA 6.317 



are present in some form in p + p data. It remains to be 
understood whether the enhancement of the ratio is due 
to parton density (multiplicity) or to collision energy. 



4 Comparison to next-to-leading order pQCD 

In this final section we discuss the improvements which 
have recently been made by next-to-leading order calcula- 
tions using more precisely parameterized fragmentation 
functions. Fragmentation functions for separated quark 
flavors have been notoriously difficult to obtain due to the 
lack of sufficiently precise collider data. However, OPAL 
has recently published flavor tagged data from e+ + e+ 
collisions which allowed theorists to compute better frag- 
mentation functions [16] . 

In figures [4] and [5] we compare two different NLO cal- 
culations to our Kg and A data . The first one (black 
lines) uses older FF by Kniehl et al. (KKP) and Vogel- 
sang et al ( WV) [17] . The second one (red lines) was done 
by Albino et al. (AKK) using more recent FF based on 
the light flavor tagged OPAL data [18] . 

Clearly, these newer parameterizations improve the de- 
scription of our A data greatly. However, in order to 
achieve this agreement, they fix the initial gluon to A 
fragmentation function (Dg) to that of the proton, then 
estimate that an additional scaling factor of 3 is neces- 
sary to achieve agreement with STAR data. However, this 
modified FF for Dg also works well in describing the p+p 
SPS data at y^s — 630GeV. It therefore appears that the 
STAR data is a better constraint for the high z part of 
the gluon fragmentation function than the OPAL e + + e~ 
data. Similar conclusions with respect to the important 
role of p + p collisions have been drawn elsewhere [19] . 




P T [GeV/c] 



Fig. 4. pr spectra for Kg at midrapity (| y |< 0.5) from p + p 
at y/s = 200GeV compared to two different NLO calculations. 
Dashed lines indicate the scale uncertainty of the NLO calcu- 
lation, ie. fi = 0.5pt (lower), p, — 2pr (upper). 
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Fig. 5. pt spectra for A at midrapity (| y |< 0.5) from p + p 
at ^J~s = 200GeV compared to two different NLO calculations. 

5 Summary 

We have shown that the theoretical description of identi- 
fied strange particles in p + p and p + p collisions is still 
not fully understood. This is especially important since 
these models are now extensively used to predict observ- 
ables for the LHC-era, and therefore one should be aware 
of their limitations. Phenomenological LO models can be 
tuned to describe the data but still struggle to describe 
baryon production at intermediate px- 

Baryon production, and in particular the baryon to 
meson yield ratios at intermediate pr, are one of the "hot" 
topics in current heavy- ion research at RHIC. The p + p 
data presented here allows us to look at the ratio in ele- 
mentary collisions and check how well it is understood in 
a simple system. The fact that PYTHIA baryon produc- 
tion parameters need to be tuned quite considerably to 
achieve an agreement is interesting. We also showed that 
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at the higher energies, i.e. y/s — 630 GeV, this difference 
is even larger. This is an indication that the baryon to 
meson effects previously observed in heavy-ion collisions 
are present in some form hi p + p data, and that the asso- 
ciated physics phenomena therefore need to be explained 
without requiring the presence of a quark-gluon plasma. 

Next-to-Leading order calculations have greatly im- 
proved with light flavor tagged fragmentation functions. 
However the high-z range of the gluon FF previously ex- 
tracted from e + + e~ data seems inconsistent with p + p 
and p + p data, indicating that RHIC data could be valu- 
able in constraining the gluon FF. 
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